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1 Introduction 



Experiments involving kaons or kaon decays have reached a high level of accuracy and are 
therefore very sensitive to probe new physics beyond the minimal Standard Model (SM). 
For example, studying the small mass difference AMk between the long-lived kaon Kl 
and the short-lived one Ks is one such measurement of amazing precision. In the K^K^ 
system, the observable AM^/Mk is very accurately established at the level of lO^"*^^ in 
line with predictions in the SM 0. Other experiments sensitive to new-physics scenarios 
look for possible lepton-fiavour-violating decays of the type Kl e/i. The experimental 
upper bound on the branching ratio B{Kl e/i) is very stringent [jll, i.e., 

5cxp (i^L ^ e/i) < 3.3x10-11, (1.1) 

at 90% confidence level. Ongoing experiments at BNL, KEK and, more recently, at DA$NE 
are expected to lower the above bound by one order of magnitude. Although the decay 
Kl efi is strictly forbidden in the SM due to the absolute conservation of the individual 
leptonic quantum numbers, lepton-fiavour violation occurs almost inevitably in many of 
its extensions. For example, many studies have been devoted to analyze such a decay in 
models based on technicolour supersymmetry horizontal symmetry ||^, extended 
scalar and gauge sector |17|,|8|J9|JT0[| . Thus, experiments on Kl — > e/i provide a unique 
opportunity either to discover new physics or to constrain efficiently the parameter space 
of these theories. 

In this paper, we shall present a careful analysis of Kl e/i in left-right models, which 
are realized by the gauge group SU(2)/i;(8)SU(2)/^ ®U(1){_b-l) and include heavy Majorana 
neutrinos [^. These theories also allow for scenarios with high Dirac mass terms of order of 
top quark mass without invalidating experimental limits on light neutrino masses. In these 
scenarios, the usual see-saw suppression relations [0 for the light-heavy neutrino mixings 



may be avoided. In particular, it has been found that the large Dirac masses give rise 



to a non-decoupling behaviour of the heavy neutrinos in loops [13,n|, thereby enhancing 
significantly the decay rate for possible lepton-fiavour-violating processes, such as Z — >■ er 
||T5| and the size of related new-physics observables at LEPl [|I6|,|T^. This feature makes 
the left-right models very attractive. Apart from taking properly into consideration the 
non-decoupling effect of heavy neutrinos, we shall also analyze the impact of other low- 
energy data on the theoretical prediction for B{Kl e/i), such as those coming from the 
non-observation of the decay /x — e7 and the KlKs mass difference. 

Another important issue, which is to be discussed in detail, is the lack of gauge 
invariance, when only contributions from one-loop box diagrams to the decay Kl —>■ efi are 



2 



considered. Studies on the K^K^ system revealed that box diagrams by themselves 

do not form a gauge-invariant set of loop graphs in the left-right models. In addition to the 
box graphs, one has to take into account self-energy and vertex corrections involving very 
heavy Higgs bosons, which mediate flavour- changing neutral currents (FCNC). Even if the 
FCNC Higgs bosons were taken to have all infinite masses, their effect would not decouple 
from the loop. In the infinite mass limit for the heavy Higgs scalars, the residual terms 
are also gauge dependent and cancel against similar gauge-dependent terms coming from 
the box graphs, leading eventually to a gauge-invariant result. While the inclusion of the 
gauge-dependent Higgs self-energies to the KlKs mass difference leads to little change in 
the value predicted by the box graphs, the situation for the decay e/i is strikingly 

different. We find that the residual gauge-dependent terms originating from the self-energy 
and vertex graphs may be comparable to or even larger than the box corrections. As 
a consequence, the constraints on the parameter space of the left-right models will now 
become much more severe than those obtained in earlier articles [§,0]. Our calculation of 
the gauge-independent decay amplitude for — > efi at one loop and the analysis of its 
immediate phenomenological consequence constitute novel aspects for the left-right models, 
which have not been studied in detail before. 

The paper is organized as follows: In Section 2, we describe the basic low-energy 
structure of the left-right models and review many useful relations among mixing param- 
eters and neutrino masses, as well as relations involving gauge- and Higgs-boson masses. 
In Section 3, we calculate the matrix element for the decay —>■ efi and derive mass 
limits on the FCNC Higgs masses from the tree-level contribution and other theoretical 
considerations. All relevant Feynman rules and one-loop analytic expressions pertaining 
to Ki — > e/i are given in Appendices A and B, respectively. In Section 4, we show how 
gauge independence gets restored, when Higgs-mediated self-energy and vertex graphs are 
taken into account in addition to the box diagrams within the well-defined on-shell skele- 
ton renormalization scheme. A detailed discussion on the dependence of B{Kl — > e/i) on 
various kinematic parameters follow in Section 5. Numerical estimates for the actual size of 
B{Ki —>■ e/i) are given within manifest as well as non- manifest left-right symmetric models. 
This enables one to deduce combined constraints on the parameters of these models. Our 
conclusions are summarized in Section 6. 



2 SU(2)i?(g)SU(2)2,(g)U(l)(5_^) model 

The idea of having spontaneous breakdown of both gauge and discrete symmetries can natu- 
rally be realized in left-right theories based on the gauge group SU(2)i?®SU(2)L(8)U(l)(B-L) 
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20| , |2T| , p!T| , where B and L denote baryon and lepton numbers, respectively. Such a real- 



ization may arise from the low-energy limit of SO(IO) grand unified theories. Specifically, 
one of the most appealing breaking patterns of SO(IO) down to the SM is the following: 

SO(IO) ^ SU(4)ps®SU(2)^®SU(2)^ -> SU(3)c ® SU(2)k ® SU(2)i ® U(1)(b_l) 
^ SU(3),® SU(2)i®U(l)y. (2.1) 

Among the several Higgs scalar representations which can lead to the SM, we shall adopt the 



field content of the left-right model introduced first in |jTT[, which includes heavy Majorana 
neutrinos. 

We will now describe the basic low-energy structure of the left-right models. In a 
three-generation left-right model, the left-handed and right-handed chiral fields are grouped 
in separate weak iso-doublets as follows: 

Li= I 1 : (0,1/2,-1) ^'«= I ) ■ (1/2,0,-1), (2.2) 




L 



u' 



g'^=| ,J : (0,1/2,1/3) ^«=|^rf/j ^(1/2,0,1/3), (2.3) 



where the weak eigenstates are denoted by a prime. In Eq. (|2.2|) , we have also 
indicated the assignment of quantum numbers to fermions under the gauge group 
SU(2)r(S)SU(2)/^(8)U(1)(b_l). Evidently, the above left-right symmetric extension of the 
SM gauge sector requires the presence of right-handed neutrinos. 

The Higgs fields introduced into the model to break left-right gauge symmetry down 
to U{l)em Q are: 



(2.4) 

where $ is a (1/2*, 1/2, 0) Higgs bidoublet, and and are two complex Higgs triplets 
with quantum numbers (0, 1, 2) and (1, 0, 2), respectively. To avoid excessive complication, 
we will assume that only (p^ and 6^ acquire vacuum expectation values (VEV's): (05) = 
and (5^) = vr/ V2. Such a scheme is also favoured for phenomenological reasons 
|2^ . At the tree level |]23|,|8[|, the vanishing of the VEV's for the neutral fields 02 and 6^ 
can be achieved by imposing invariance of the left-right symmetric Higgs potential under 
the discrete symmetry D: $ ^ z$, —>■ Al and A^ — A/j. As has been observed 
in p4|, however, the D symmetry may suppress unnaturally large FCNC's but leads to 



zero masses for the down quarks and the charged leptons. Even though the latter may be 
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viewed as a viable approximation for Kl — > e/x, the presence of non-zero m^^ and m/^ would 
require infinite renormalization of the VEV of 5°. To overcome this problem, the authors 
in suggested that stability of such a model under high order quantum corrections can 



naturally be achieved within a D-parity breaking scenario. Since we shall treat charged 
leptons and down-type quarks as being strictly massless at the quantum level, the above 
problem does not occur in our calculations. 

The left- right scenario of our interest is further described in Ref. |2^, under case 
(d). The fact that (6^) = {(f)^ = in this case implies the absence of mixing between the 
charged gauge bosons W^^ and Wp^. Thus, W^^ may be identified as the known boson 
with mass My/ observed experimentally, whereas is an extra gauge boson having a mass 
Mr larger than at least 0.5 TeV [0]. In addition, the minimal left-right model predicts two 
massive neutral gauge bosons, Zl and Zr, which have a non- vanishing tree-level mixing of 
order i^\/v\ ~ 10~^. Nevertheless, Z/, /^-mediated interactions are not of much relevance 
here, since they do not directly enter our discussion of ep,. 

The scalar spectrum of the left-right model consists of 20 degrees of freedom, of which 
six degrees become the longitudinal components of the gauge bosons: W^, and 
Zji. From the remaining 14 physical states, only the charged scalars h^, the CP-even Higgs 
boson 02 = ^e(f)2/V^ and the CP-odd Higgs scalar 02 = Q'm^^/v^ are of direct relevance 
to the decay Kl — > e/i. The field is the orthogonal component of the right-handed 
would-be Goldstone boson G^. The fields and are not pure mass-eigenstates but 
admixtures of (pf and 6^, viz. 

As we will see in Section 3, the scalars 02* must be very heavy, i.e., heavier than 10 TeV, in 
order to avoid large FCNC contributions to the KlKs mass difference. Furthermore, after 
diagonalization of the relevant mass matrices, /i^ and 02 * are found to be nearly degenerate 
H. In fact, one obtains 



Ml, = Ml, + 0{M^^) = Ml + 0{M'^) (2.6) 



''2 



Therefore, it is reasonable to assume that the Higgs scalars h , 02 and 02 are exactly 
degenerate. 

The left-right model admits the presence of S — L-violating operators in the Yukawa 
sector, which are introduced by the triplet fields ^l,r in the following way: 



C^^ = -^^|^(/.,,l£.,,A,L',^ + Lg..,A«L'^J + H.c, (2.7) 
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where Eij is the usual Levi-Civita tensor, mMij are Majorana mass terms and hij are Yukawa 
couphngs, which are much smaller than unity owing to phenomenological constraints from 
muon and r decays |^ . After spontaneous symmetry breaking, the Lagrangian describing 
the neutrino mass matrix is given by 

^Lss = -I (^L ^'^) ( ) + H-c. 

= ulM^hr + H.c, (2.8) 
where the 6 x 6-dimensional matrix takes on the known see-saw-type form, i.e., 

M- = ( \ ] . (2.9) 

y m|) niM J 

The Dirac mass matrix m/) is a 3 x 3-dimensional matrix originating from the breaking of 
the SU(2)i sector. In general, the matrix M" is always diagonalizable by a 6 x 6 unitary 
matrix through 

W^M'^W = M\ (2.10) 

After diagonalization, one gets six physical Majorana neutrinos n^. Three of the six mass- 
eigenstates are assigned to the ordinary light neutrinos, u^, and z/^, while the remaining 
three states, A^i, N2 and A^3, describe new Majorana neutrinos, which must be heavier than 
about 100 GeV for phenomenological reasons [ p6[| . 

By analogy, the 3x3 charged-lepton mass matrix can always be diagonalized 
through the bi-unitary transformation 

VlM^VI = M', (2.11) 

where Vl and Vr are unitary matrices that relate the weak fields and Vr to the mass 
eigenstates li and Ir, respectively. Adopting the conventions of Ref. [^, we can now write 



the charged- current interaction of the Majorana neutrinos and charged leptons li with 
the left-handed gauge boson (= W^). The charged-current Lagrangian is given by 

CZ' = -^W^^Bth.^LU, + H.c, (2.12) 

where is the weak coupling constant, = 1 — = M^^/M^, Pl(r) = [1 ~ (+)75]/2 
and is the 3x6 matrix 

fc=i 
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which describes the mixing of the left-handed charged leptons with the neutrinos. Further- 
more, the Zi^ boson couples to the Majorana neutrinos through the 6x6 mixing matrix 

c% = Y^ ui:m;- (2.14) 

k=l 



In an analogous way, the mixing matrices -B^f and Cf^ 



I] 1 

BS =EV^i?.f/r, and = Y: [/,7f/i;, (2.15) 

fc=4 fc=4 

mediate the interactions of the and Zr bosons to the charged leptons and Majorana 
neutrinos. In the left-right models, the quark sector is modified as well, namely, there 
are additional couplings of the fermionic fields to the gauge and extra Higgs bosons. All 
Feynman rules that are necessary for the calculation of — > e/i are relegated to Appendix 
A. 

The left-handed flavour- mixing matrices and C^, which are identical to the re- 
spective mixing matrices B and C of the SM with right-handed neutrinos, satisfy a number 
of identities that may be found in p6|JI^ . In addition, the right-handed mixing matrices 
B^ and obey the relations ||17||: 



Ei^^.i?.^, = 0, E^fSj = Cj, C^* + C5 = 6,,. (2.16) 
1=1 1=1 

These identities is a direct consequence of the unitarity of the theory, assuring its renor- 
malizability. To simplify further the analysis of e/i, we shall assume that only the 

electron and muon families mix effectively, whereas the r family stays unmixed. In such 
an effective two-generation-mixing model, one obtains ||T3| , p!^ 

Bk = ' ^k = -y^^ ^ (2-17) 



where the ratio p = mj^^/mj^^ is always understood to be greater or equal than unity and 
s'l (with I = e,fi) are the light-heavy neutrino mixings introduced in i.e., 



ElBmf ^ [^D-r^D] . (2.18) 



i=l \ 

In addition, for the mixings C^.^., one flnds 



^N-iNo — — "^P^^^ ^NoNo ■ (2.19) 
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Finally, with the help of Eq. ( p.l6| ) together with Eq. ( p.l5| ), we can find the mixings B 



up to an arbitrary Cabbibo-type angle 9^ [|T7|. The mixings Bjj^, are given by 



nR 



TjR 



COS On 
— COS Or 



sin 9 



R 



1 — 

^N2Nx 



nR 



— i sin 9 



R 



C-NiNi)i^ CN2N2) I^'aTjAtJ 



1 



i cos 9r 



1 ^NiNi 



1/2 



\^NiN2 I 



1/2 



.(2.20) 



Counting now the number of independent kinematic parameters present in the leptonic 
sector of our two-generation left-right model, we find five free quantities: the masses of the 
two heavy Majorana neutrinos, m]\r-^ and 171^^ [or equivalently miy-^ and p], the two light- 
heavy neutrino mixings, (s^")^ and (s^^)^, which are severely constrained by low-energy 
data as we will see later on in this section, and the angle 9r. 

It is now interesting to remark that M'^ in Eq. (p.9|) reduces to the traditional seesaw 



scheme |jT2| in the limit niM ^ "nrio- In the see-saw mass pattern, the light-heavy neutrino 
mixings s'l are very suppressed, since they are governed from the relation s'l ~ ^Jm^Jmi^. 
The usual seesaw models predict very small rates for lepton-fiavour-violating decays, such 
as Ki — i> e/i, /i — i> eee, etc., which are beyond the realm of detection in any foreseeable 
experiment. The situation may change drastically if a non-trivial mixing between two 
families ||2^ is introduced. To give an example, consider the following Dirac and Majorana 
sub- matrices in Eq. ( p.9|) : 



mo 




rriM 




(2.21) 



where a, b, A and fi are arbitrary complex numbers. It is then easy to verify that the rank 
of the 4x4 matrix M'^ is two. This implies that two eigenstates are exactly massless at 
the tree level, whereas the other two are massive. The scheme in Eq. ( |2.21| ) may also be 
derived from an horizontal symmetry, which is broken softly by a lepton-number-violating 
operator proportional to p. In fact, if /i arises from the spontaneous breakdown of a global 
symmetry, then it should be p -C A in order to avoid large Majoron couplings to electron 
and u, d quarks, as is required from astrophysical considerations . In this scenario, the 
light neutrinos acquire masses radiatively in agreement with experimental upper bounds 
p6| . Most importantly, the mixings s^l now scale as ~ a/ A, b/A and therefore are large. 
In particular, we should note that the derivation of the mixings B^ and B^ in Eqs. ( |2.17|) 
and (|2.2CI|) are based on non-seesaw scenarios of the generic form (|2.21|) . These mixings 
can only be constrained from a global analysis of low-energy data. 
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There are many experimental data that can be used to constrain the hght-heavy 



neutrino mixings |27|. Here, we shall however focus only on those limits that involve the 
e and fi families and are relevant for the decay Kl —>■ efi. We shall also assume that 
the two heavy Majorana neutrinos A^^i and N2 are almost degenerate in mass, namely, 
rriNi = = ^N- Thus, charged current universality in tt decays and in the decay 
fi evv lead to the limits: 

{sl?^ {s7f < 0.010. (2.22) 

Furthermore, in the large mjv limit, heavy Majorana neutrinos give rise to fi —>■ decays 
with branching ratio 

5(/i-e7) = ^-^{sl^nsTf. (2.23) 
On the experimental side, Bcxp{fi 67) < 4.9 x 10^^^, which yields the tight constraint: 

(sl)'(sI")' < 5.6 X 10-^ (2.24) 



For the model at hand, the absence of /x — e conversion events in nuclei |Q may lead to 
more stringent bounds due to the non-decoupling behaviour of the heavy neutrinos in the 



Zefi coupling |T^. A first estimate gives 

Iststlistr + isln'] - 10'^ X ^. (2.25) 



N 



If the heavy neutrinos have TeV mass, both constraints (|2.24|) and (|2.25|) are then com- 



parable. For very heavy neutrinos of order 10 TeV, the bound due to /i — e conversion is 
more restrictive. Of course, one can completely avoid the latter constraint by assuming, for 
example, that the muon number is practically conserved by the left-handed interactions, 
which merely means that 6 <^ a in Eq. ( |2.21| ). Even if one takes = 0, Kl — > e/i can 
still proceed via the exchange of and bosons in the loop. Finally, we find that the 
decay fi —>■ eee provides in general a weaker constraint than that given in Eq. ( p. 251) . 



In the see-saw type matrix ( p.9| ), the Dirac and Majorana mass terms, m^j and mM, 
cannot be arbitrarily large but are constrained from above by triviality and perturbative 
unitarity bounds. Renormalization-group-triviality bounds, which are mainly controlled 
by large contributions to the quartic scalar couplings in the Higgs potential, usually lead 
to the limit: uid < 2Mw / Qw ~ 300 GeV. However, the coexistence of large Majorana 
masses tum ^ m^) in the loop of the Higgs potential will screen the one- loop contributions 
of the heavy neutrinos to the Higgs potential considerably. In fact, the low-energy data 
mentioned above require mo/TaM < 10^^. Therefore, a more reliable constraint arises due 
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to perturbative unitarity, where the inequahty r(A^ ^ ^W^l", ^Zl) < 111^/2 is imposed. 
This inequahty amounts to the bound 



(2.26) 



Similarly, the Majorana mass itim may be constrained from an analogous inequality r(A^ 
IW^, vZr) < mAr/2, which gives the bound 



< 



niM ~ 10 Mr. 

The limits derived in Eqs. (^.26|) and (|2.271) are used throughout our analysis. 



(2.27) 



3 Matrix element of Kl e/i 

In the left-right models, the lepton-flavour violating decay Kl —>■ efi can occur at the 
tree level via the exchange of FCNC Higgs scalars, as shown in Fig. 1. However, for 
very heavy FCNC Higgs masses, the tree-level contribution gets rather suppressed due to 
the intermediate Higgs-boson propagator, while one-loop graphs shown in Figs. 2-4 can, 
in principle, become very significant, since they do not depend explicitly on the FCNC 
Higgs masses. Therefore, we shall include in our study all the relevant one-loop diagrams 
contributing to Kl — > e/i. 

The one-loop matrix element of the decay Kl — ^ e/i has the general form 

r(i?°^eV) = (i^(^) j^mi.PLs\K')u,rPLVe{A + vB), (3.1) 

where the reduced amplitudes A and B are given in Appendix B. The parameter rj in Eq. 
( p.l|) is an enhancement factor originating from chirality-flipping operators which enter in 
the kaon-to- vacuum matrix element. The factor rj is given by 

{0\drrPRs\K'')u^-f^-f,PLVe 



V 



{0\drPLs\K^)u^.laPLVe 



^ t ^ ^50, (3.2) 

where is the mass of K^, and and denote the current d- and s-quark masses, 
respectively ||3T[. In the derivation of the parameter rj in Eq. ( |3.2| ), we have taken into 



account the PCAC hypothesis and the fact that 

{0\d^,luPRs\K\p)) = g,,{0\d^,s\K'^{p)) + ^(0|rf[7^,7.]753|ir°(p)). (3.3) 
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The second term on the RHS of Eq. ( |3.3| ) must be a Lorentz tensor antisymmetric in p^p" 
and hence vanishes identically. Applying the PCAC hypothesis to the hadronic matrix 
element 

{0\d{x)^,^,s{x)\K\p)) = ip,fKe-'^\ (3.4) 
where fx is the kaon decay constant, we can calculate the matrix element 

(0|J(x)75^(x)|ir°(p)) = -I e-^^-. (3.5) 

Using the relations ( p.3|) , (|3.4|) , ( ^75|) and little algebra, the second equality in Eq. (|3.2|) 
follows easily. 

Substituting Eq. ( |3.4|) into Eq. ( |3.1| ), we find the branching ratio 

B{Kl ^ e/i) = 4.1 X 10"^ 1^ + r/-B|^ (3.6) 

where the charge of the leptons in the final state is not tagged. One arrives at the same 
result for B{Kl — * e/x) in Eq. (|3.6|) , by using relations based on isospin invariance between 
the decay amplitudes of —>■ /i~e+ and f^^t^a- 



^2 



Fig. 1: Tree-level graph contributing to — > 

The tree-level (j) 2^ -exchange shown in Fig. 1 gives rise to a contribution to T{K^ 
efi), which is proportional to the reduced amphtude 



r 1 



(3.7) 



whereas Atree = 0. In Eq. (p.7|), we have defined Ac,t = m^^/M^, Aq = m^^/M^, Xji 
MlrjMl and A, = Mj./M^,. 
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— * — rwvvi — *— 



— * — rvwvi — *- 



d h+ 



Ui 



(a) 



(b) 



'"Tin 



h~ 
(c) 



Fig. 2: Gauge-independent sub-set of Feynman graphs (group A) contributing to Kl e^u 




(g) (h) 

Fig. 3: Gauge-independent sub-set of Feynman graphs (group B) contributing to Kl efi 



There are tight experimental constraints on the masses of the FCNC scalars and the 
mixing angles, which depend on the flavour structure of the left-right model. For example, 
in non- manifest left-right models 0, the experimental limit coming from the mass difference 
between Kl and Ks leads to the lower mass bound 

M^r ~ M^. ~ 9 TeV. (3.8) 

If we impose a manifest left right-symmetry on the model, i.e., = = V, the mass 



bound (p.8|) is even more severe, and the FCNC scalars must be heavier than 30 TeV |32 
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On the other hand, in the manifest left-right symmetric model, the tree-level contribution 
to Kl e/i yields the following constraint: 

M > r ( ^ ( \ ( ^ ( ^ TpV 

"^^J'^'H^J iaoij lisoGivj lioTwj ^'-'^ 

which is comparable to the one obtained from the KlKs mass difference. 

Another important bound comes from the requirement that the validity of perturba- 
tive expansion be preserved. The unitarity bound derived in this way is translated into the 



approximate inequality ||33 



M^r^r ~ 15 Mr. (3.10) 

The above bound will be considered throughout our phenomenological analysis in Section 
5. 




(g) (h) 

Fig. 4: Gauge-independent sub-set of Feynman graphs (group C) contributing to Kl — > 



If the FCNC Higgs scalars are sufficiently heavy, e.g. much heavier than 10 TeV 
according to Eqs. ( ^.8| ) and (^.9|), the tree- level contribution to T{K^ e/i) may become 
smaller than the present experimental limits. However, the one-loop box diagrams shown 
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in Figs. 2-5, being independent of the heavy neutral Higgs masses, can still give significant 
contributions to the decay efi. Therefore, the impact of the loop corrections on 

the phenomenological predictions should be studied very carefully. Thereby it should be 
kept in mind that the box diagrams alone do not form a gauge-invariant set, but vertex 
and Higgs vacuum polarization graphs must be included. This issue together with that of 
renormalization will be discussed in detail in the next section. 



d 



R 



d h+ 



l^/\/V\>r 



Ui 



rir 



Ui 



rir 



s h- H~ 

(a) 



— > — ^AAAJ — >— 

s /i^ 

(b) 



Fig. 5: Gauge-independent sub-set of Feynman graphs (group D) contributing to Kl e/u 



4 Gauge independence and renormalization 

In general, the and propagators depend on the gauge used to remove the unphys- 
ical degrees of freedom and on the parameter chosen for fixing such a gauge. There are 
many gauges to carry out this procedure, such as covariant gauges, non-covariant axial 
gauges, etc. Hence, the and propagators may introduce gauge-fixing-parameter 
dependence if only an arbitrary sub-set of loop graphs is considered. For example, con- 
sidering box graphs only to describe the K^K^ mixing in left-right models is found to be 
insufficient to cancel such a dependence on the gauge- fixing parameter ^ [^|TH|,|T^3. 

In the following, we shall closer examine the gauge dependence of T{K^ — > 
on the Wl and Wr propagators in the covariant gauge. The R^ gauge may be viewed 
as the most practicable class of gauges endowed with the manifest properties of Lorentz 
invariance and renormalizability. In particular, we wish to identify all those diagrams, 
which, together with the box graphs, form a minimal gauge-independent set. To this end, 
we have divided all the radiative corrections into four groups A-D, depending on the way 
that the Wr and Wl propagators enter in the loop. These assignments are represented 
in Figs. 2-5. It is then not difficult to verify that the groups of box graphs A and D are 
separately independent of the gauge-fixing parameter ^. In the other classes of graphs, i.e., 
B and C, however, one must consider additional diagrams containing FCNC Higgs scalars 
and all relevant tadpoles, which are not shown in Figs. 3 and 4. Since the self-energy 
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diagrams in groups B and C are UV infinite, we will simplify our task by proving gauge 
independence after renormalization. This is indeed a rather economical way to check ^ 
independence, since a vast number of tadpoles, being momentum independent, will drop 
out, when all relevant renormalization subtractions are performed. For that purpose, we 
shall adopt the on-shell skeleton (OSS) renormalization scheme, which is also discussed in 
Ref. m. 



The main virtue of the OSS scheme is its explicit maintenance of gauge independence 
during the process of renormalization. The OSS renormalization is based on subtract- 
ing lower n-point correlation functions from the transition amplitude under consideration. 
The lower n-point correlation functions, which normally represent self-energy and vertex 
graphs, are formally gauge independent, since they are evaluated at on-shell external mo- 
menta. More explicitly, the OSS scheme may be described as follows. Whenever we have 
to renormalize a vertex sub-graph, e.g., the one-loop 02/^~e"'' coupling r0r(p^) shown in 
Figs. 4(c) and 4(d), we have just to make the subtraction 

rS(/) = r^5(p') -r^5(M|). (4.i) 

The above operation is called Rl subtraction. By analogy, we renormalize all the self- 
energy graphs involving the FCNC Higgs bosons, e.g., n^r(p2^ Figs. 4(g) and 4(h), by 
making two subtractions defined as 

nf (/) = n,.(/) - n,.(M4.) - (p^ - m|) An,.(p2) , (4.2) 



■^2 



which is characterized as R2 operation. In our analysis, we shall not study possible effects 
due to renormalization scheme dependence [|4], which are of high order. In any case, the 



lepton-flavour-violating decay Kl — > e/x would not constitute the best place to look for 
such effects. 

In the following, we shall see how the box graphs in class C (Fig. 4) are ^ dependent 
and how gauge independence gets restored if the relevant Higgs self-energies and vertex 
graphs are taken into account in the OSS renormalization scheme. The same procedure 
may be applied to the class B. The class A and D are separately gauge independent and 
UV finite. Obviously, OSS renormalization cannot apply to the groups A and D, since they 
only consist of box diagrams. 
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Fig. 6: Diagrammatic definitions of tree-level amplitudes pertinent to class C. 



4.1 Gauge dependence of box diagrams 

We start examining the gauge dependence of the box diagrams in class C shown in Fig. 
4 in the gauges. In this class of gauges, the gauge-boson propagator, e.g., of the Wl 
boson may conveniently be decomposed as follows: 

A,.(^) = U,.{k) - ^A,{k) , (4.3) 

where 

is the respective VFt-boson propagator in the unitary gauge and 

1 

is the ^-dependent part of the VF^-boson propagator. Note that A^(/c) coincides with the 
propagators of the would-be Goldstone boson Gl and the ghost fields c^, c^, which cor- 
respond to Wl- Similarly, one can perform an analogous decomposition for l^^-boson 
propagator, which will depend on another gauge-fixing parameter, e.g., Such a decom- 
position is very useful to study the gauge independence of the class B, which is, however, 
conceptually very similar to that of the class C considered here. Finally, the propagator 
for the free charged Higgs bosons h"^ is given by 

D.ik) = ^ , (4.6) 
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which is independent of ^. 

To make use of the propagator decomposition mentioned above, we write the amph- 
tude for the diagram (a) in Fig. 4 as a sum of two sub-amphtudes: A^a = + -M^- 
Omitting integration over the loop momentum, the sub-amphtudes are given by 

= -T^\k)fl:'^{-k)A^^{k)Dr,{k') , (4.7) 
j^G ^ -T^h^k)f^'^{-k)Ai:{k)Dh{k') . (4.8) 

where k and k' are the loop momenta carried by the and charged Higgs bosons h~ , 
respectively. In Eqs. (^4.7| ) and ( [4 .81) , T^'*, T^^, and T^'-' are the tree-level amplitudes 
depicted in Fig. 6. These tree- level amplitudes read: 



1 



2 



r. ^,B^:FLVe{p[) , (4.9) 



+ /if - 



2 



T^\p[,p'„k) = -^(^j MP2)[c,m,BXPR+[sl5,a-Ci^)m,B;^^F, 

¥ B^:{^ePR - m^FL)veip[) , (4.11) 

x(m,Pfi-m,Pi)V^,fM,(pi). (4.12) 
Using the propagator decomposition of Eq. (|4.3|) , one can now write 

= Ml^Mi, (4.13) 

with 

-^a = -T^\k)f!:'^{-k)U^''{k)D,,{k') , (4.14) 

Mi = -^T;^\k)f!:'^{-k)A^{k)D,{k'). (4.15) 

Here, Ai^ denotes the amplitude of box diagram 4(a) in the unitary gauge, which does 
not explicitly depend on the gauge parameter ^. In the unitary gauge, the unphysical 
Goldstone bosons (and ghosts) are absent. 
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The tree-level amplitudes given in Fig. 6 satisfy the following Ward identities: 

" -Tj^'^ik) = -T^'^-El, (4.16) 



Mw 

-^T;^(A;) = T^'^ + Er, (4.17) 
-^T;^(A;) = T^^ + En, (4.18) 



iw 

f^^{k) = -T^'^-^L, (4.19) 



Ml 



where we have defined 

^L(pi,P2) 



Er{pi,P2) = M^l 



\ 2 



\ 2 



C/3 


^^(p'2)i?« 5i:m„Pi^;e(p;) , 


(4.20) 




Vdip2)V^/VimiPRUs{pi) , 


(4.21) 


Ml 




(4.22) 






(4.23) 



El{pi,P2) = 

With the help of the above Ward identities, the amphtude Ml can then be rewritten as 

Mi = -M': + B^ + B^s, (4-24) 

where 

= A^{k)Dh{k')(T^''ER + ELf""^) , (4.25) 

S§ = A^{k)Dh{k')ELER . (4.26) 

Following a similar procedure for the diagram 4(b), we obtain 

Ml = -Mf + B^y + Bl. (4.27) 

The final expression for the box graphs 4(a) and 4(b) takes the form 

Albox = Ma + Mb 

= M'^ + M^ + Mi + Ml + M^ + Mf 

= Ml^ + Bv + Bs, (4.28) 

where M^^^ — M^ + M^ is the contribution of the box diagrams evaluated in the unitary 

gauge, and 

Bv = A^{k)Dh{k') (t^'^Er + ElT^^ + T^'^El + ErT^^) , (4.29) 

^5 = A^{k)Dh{k')(ELER + EREL) . (4.30) 
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The presence of the terms By and Bs, which are proportional to the gauge- dependent 
propagator A^{k), imphes that the box diagrams in class C do not form a gauge-invariant 
set, and that additional vertex and self-energy graphs must be considered. 



4.2 OSS renormalization of vertex diagrams 

We now consider the contribution of the vertex diagrams 4(c)-(f) in the OSS renormal- 
ization scheme. After performing the Rl subtraction given in Eq. ( [4.1|) , we obtain the 
following matrix element for the vertex graphs: 

^=^r,i P ^^^^ P ^^^ip 

where we have defined 

rj.(p;,p'2) = - ^Es{p\.p',)-, ^\{P\.P'2) = -'^Ep{p[,p',), (4.32) 

^ 9wCi3 ^2 g^Cj3 

rUp^,p,) = -M}Les{p^,p2); rUp^,p,) = -'-Mkep(p„p,) , (4.33) 

with 

EsAp'i^P2) = En{p[,p'2)±EL{p[,p',) , (4.34) 

Es,piPi,P2) = Er{pi,P2) ± El{pi,P2) ■ (4.35) 



The functions El, Eji, Ept, and El have previously been defined in Eqs. (|4.2CI|) - (^4.23|) , 
respectively. 

We now proceed by decomposing A^^^tex terms of an _Rl-subtracted amplitude in 
the unitary gauge and the remainder, i.e., 

-^vertex = -^vertex + Vs + Vb . (4.36) 

In analogy with Eqs. ( |4.7| ) and ( |4.8| ), the one-loop vertices involving the virtual states 
Wl and Gl may also be separated into two terms: T^, = + and for T^, likewise. 
Employing the identity 



Ml -Ml 



Rl 



= ^r=%r=(/) - [iL^r°(/) , (4.37) 



we can cast the different terms on the RHS of Eq. ( |4.36| ) into the form 

•'^Ivertex "T Y V ]lf2 _ /Kf2 V ^ Jlf"^ — M"^ 



h if h <fi 
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+ 



B 



p2 - M, 



Ml 



m 



M2 







Ml ^ 



I f° 



(4.38) 
(4.39) 



Using the Feynman rules given in Appendix A, we can calculate the analytic expres- 
sions for the one-loop vertices. These are given by 



rg(/) 
rg(p^) 



w 



2M, 
2Mw 

2Mw 

gwCfs 

2Mw 
2 



{Ml-Mlr)A^{k)D,{k'){T^^ + T^^ 



{Ml-Ml)A^{k)Dn{k')(T 



nhG rpGh 



{Ml-Ml.)A^{k)D,{k'){f^^ + f 
{Ml-Ml^A^{k)D,{k'){f'^^-f^^' 



(4.40) 
(4.41) 

(4.42) 
(4.43) 



{Ml-Ml)\u^'^{k) 



{Ml-MlA U^^{k) 



'mI; 

kf^k" 



A^{k)] D,{k') (t^ - Tf {k' + p)„ (4.44) 
A^ik)] Dn{k') (T^+rr) {k' + P)„ (4.45) 



..2^^......^ ^'^'-A^(A:)]D,(A:')(Tr-Tr)(A;' + p)„ (4.46) 



2 



{Ml-MlA U^^{k) 



Ml 
k^'k' 



Ml 



A^{k)] D^{k') (t'^+ff^) {k' + p),. (4.47) 



With the help of these expressions, it is now straightforward to show that 



Vb 



-B 



V ■ 



Evidently, in the evaluation of the total matrix element, Vb will cancel against 
dependent term By coming from the box diagrams (c/. ( [4.29| )). 

We will now determine the analytic form of the ^-dependent term Vs in Eq. 
For this purpose, it is useful to introduce the abbreviation 



(4.48) 
the ^- 



r.(/) - T^iP 



We can now express T^r^p"^) as follows 



2, p' - Ml a 
Ml -Ml ^ 



rS(p^) . 



(4.49) 



'^-''ip'-Ml)EsA,{k)D,{k'). 



2Mi 



w 



(4.50) 
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Here, r^r(p^) denotes the one-loop 02/^~e+ coupling, evaluated in the unitary gauge. Since 
the second and the third term vanish under the Rl operation, as they do not depend on 
p^, we then obtain for the OSS renormalized vertex function 



2M, 



w 



ip' - Mi)A^{k)D,ik') 



Rl 



(4.51) 



The first term on the RHS of Eq. ( [4.51| ) is needed to obtain A^^rtL iii Eq. ( [4.38| ), whereas 
the remainder is a function that depends on ^ explicitly. Including both 02 and (f)\ contri- 
butions, we obtain the ^-dependent function Vs 



Vs 



ip' - Ml)A^{k)D^{k')Y\EsEs , \,2 - EpEp 



p^-MlJ- 



(4.52) 



In the following, we shall show how the ^-dependent terms Vs and Bs-, defined in Eqs. (|4.52| ) 
and ( |4.3CI|) , respectively, will cancel against corresponding ^-dependent terms coming from 
the self-energy graphs, shown in Figs. 4(g) and 4(h). 



4.3 OSS renormalization of self-energy diagrams 

We shall now calculate the amplitude related to the self-energy graphs (g)-(h) in Fig. 4. 
After carrying the R2 subtraction defined in Eq. ( [4.2| ), we may conveniently write the 
self-energy-like amplitude as 

1 



A// ^2 



1 



V=<t>2 



(4.53) 



By analogy, we also write the self-energy of the FCNC scalars as a sum of two terms: 



K(p") + K^p 



(4.54) 



where 11^ (p ) and 11^ (p ) are respectively the graphs with Wl and Gl bosons in the loop. 



The action of R2 operation on the self-energy 11^ (j9^) gives 



n 



Mi - Ml 



\MI-MIY 



ip' - M2)nJ(/^ 



Rl 



+ 



\M? - 



(4.55) 



With the help of Eq. ( [4.55|) , the self-energy- like amplitude may be decomposed as follows: 

(4.56) 



-^sdf ' + Sb + Sv 
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where 



-'^'sclf 



Sb 
Sy 



- E r° 
E K 

¥'=02 



-Ml 



A/f2 

:^)nJ(/) + n-(p^ 



iJ2 ~„ 



(M2-M2)(p2_M2) 



(4.57) 
(4.58) 
(4.59) 



We should now show that A^sdf 

in Eq. ( ^.57| ) is indeed the OSS renormahzed self- 
energy-hke amphtude in the unitary gauge. To make this exphcit, we first calculate the 
analytic expressions for all the scalar self-energies 11^ (p^) and 11^ (p^) before applying 
the R2 subtraction. Using the Feynman rules listed in Appendix A and omitting loop 
integration, the different self-energy contributions read: 



ng(p^) 



24 (9 



(f) {Ml-Ml^^/\,{k)D,{k') 



24 (9w 



{fy{Mi-Mi.rA,{k)D,{k') 



kfj_ky 



A^{k) Dt,{k'){p + k'Y{p + k') 



(4.60) 
(4.61) 

(4.62) 
(4.63) 



Taking these last results into account, we obtain for the expression within the square 
brackets, which appears on the RHS of Eq. ( |4.57| ), 



Mi 



(k'^-Ml) A^ik). (4.64) 



The first term in Eq. ( [4.64| ) is the 02 self-energy in the unitary gauge, whereas the second 
term, being only linear in p^, will vanish after the R2 operation, which requires two sub- 
tractions. As a result, A^^'if^ is indeed the OSS-renormalized self-energy-like amplitude in 
the unitary gauge. 

We are then left with the gauge-dependent terms 5*^ and 5*^, which can be conve- 
niently re-expressed as follows: 



Sb = -A^{k)Du{k')(EREL + ELER) , 
Sv = [ip'-M',)A^ik)D,ik')]''\EsEs 



Ml 



-EpE 



p-^p- 



Ml 



(4.65) 
(4.66) 
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By simple comparison of Eqs. ( [4.65| ) and ( [4.66| ) with Eqs. ( [4.30| ) and ( |4.52| ), it is easy to 
see that 

Sb = -Bs, Sv = ~Vs (4.67) 

This very last result concludes our proof of gauge-invariance of the diagrams in group C, 
presented in Fig. 4. 

To sum up, we have shown that the total amplitude Ai of the diagrams shown in Fig. 
4 is independent in the OSS renormalization scheme. As it should, A4 equals the result 
obtained in the unitary gauge, i.e., 



M 



-Mbox + -Mvertex + M^l, 
^^Ibox ^ •'^Ivertex ^ •'^Is 



'self 



(4.68) 



In particular, we have seen how the gauge-dependent terms originating from the box graphs 
cancel against self-energy and the vertex contributions after OSS renormalization. A line 
of similar steps may be followed to show the gauge independence of the group B. Therefore, 



expressions analogous to Eq. ( |4.68|) are to be used in our phenomenological considerations 
in Section 5. 



5 Discussion and numerical results 

The branching ratio B{Kl e/i) given in Eq. ( |3.6|) depends on the reduced amplitudes 
A and B, which in turn depend on many kinematic parameters. Since our analytic results 
are shown to be gauge independent in Section 4, we shall henceforth adopt the Feynman- 
't Hooft gauge in the discussion that follows. First, we will qualitatively estimate the 
dominant one-loop contributions to the chirality enhanced amplitude B. Then, we will 
compare these estimates with those obtained for the K^K^ system and underline the crucial 
difference of the one-loop results for e/i with those found for the K^K^ system. In 

addition, we will point out the main improvements of our analysis, compared to earlier 
studies. Using the complete analytic expressions listed in Appendix B, we will present 
numerical predictions for B{Kl — > e/i) in manifest and non-manifest left-right models and 
determine the allowed parameter space, when the experimental limit -Bexp in Eq. ( |1 . 1|) is 
implemented. 

We shall now discuss the qualitative behaviour of the individual contributions to the 
decay amplitude T{K^ e"*'/i~) in the limit mNi,mN2 ^ Mr ^ My/- In this heavy- 
neutrino limit, the reduced amplitude. All, which describes the contribution from the two 
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Wl bosons in the loop, behaves as 

In deriving Eq. (pTT|), we have used the approximation: s^rriN ~ "mz?- We see that 
decreases in magnitude for large rriM values when m^) is kept fixed, whereas the Dirac 
masses m^i do not decouple for fixed values of mM- Precisely the existence of this non- 
decoupling window due to high Dirac masses rriD has led to enhanced predictions for lepton- 
fiavour- and universality-violating processes in minimal extensions of the SM with right- 
handed neutrinos [^,|r^,|T^,|T^ . For ultra- heavy neutrinos, the reduced amplitude Arr 
originating from two Wr bosons shows up a behaviour quite analogous to All, 

Arr ^ - V,, V,, ^J^^^^ j • (5-2) 

In analogy to All, Arr vanishes when niM — > oo. 

As has already been noticed in |jlO|, the dominant one- loop contribution to T{K^ 



e"'"/i~) comes from the box graphs with virtual Goldstone bosons and in Figs. 
3(a) and 3(b). In Section 4 we have seen however that the box graphs are not gauge 
independent by themselves, and it is therefore important to include the respective gauge- 
dependent complements originating from graphs in Figs. 3(g) and 3(h). If these additional 
contributions are included, we observe that the gauge-independent reduced amplitude r]B is 
still chirality enhanced due to the simultaneous presence of left- and right-handed currents 
in the loop. Expanding to leading order in the large parameters (3\n = mj^/M'^ (with 
(3 = M^/Ml) and A* = mf/M^, we find 

vB{3a + 3b + 3g + 3h) ^ ^ iV,^/V,fB;!j,B^^ + Vi*V,':B;:^B^*)P{XANy^' 



X 



(l + pXtlnp + 



Xt-\ 

rL*^rR , xrR*^rL\^trnp 

In the last approximate equality of Eq. ( |5.3| ), we have used again the fact that s^'m^v ~ ttlli 
and P <^ 1. In contrast to All and Arr, rjB does not vanish in the itim/vr oo. This 
is a novel consequence of the left-right models and cannot occur in other models, in which 
gauge interactions conserve chirality, such as SU(2)i(8>U(l)y scenarios with right-handed 
neutrinos . 

Apart from restoration of gauge independence, it should be stressed again that the 
one-loop self-energies of the FCNC Higgs scalars play a very crucial role in the phenomenol- 
ogy of Kl — > e/i. In the Feynman-'t Hooft gauge, the dominant contribution to the reduced 
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amplitude B is due the Higgs self-energies in Figs. 3(g) and 3(h). Their contribution relative 
to the box diagrams may be determined from the ratio 

B{3a + 3b) (3\t ln/3 + 81n/j/Ajv 
r-box/FCNC = g(3^ + 3/,) ~ -l + lln(M2/M|) ' ^^^^ 

We see that the contribution of the FCNC Higgs self-energies to T{K^ —>■ e~^fi~) is con- 
siderably enhanced by a factor ~ 1/(3 in comparison with that of the box diagrams. From 
Eq. ( p.4|) , we can readily deduce that rbox/FCNC < 1, unless P < 0.01 and < 1 TeV. 
However, for the parameter range, in which the box graphs are dominant, we find that the 
branching ratio is at least by one order of magnitude smaller than the present experimental 
limit -Bexp given in Eq. (|1.1| ). 

In the manifest left-right symmetric model (LRSM) {i.e., = = V), the chirality 



enhanced amplitude rjE in Eq. ([5.3[) can lead to a combined constraint on m^) and M^, 



when the experimental limit on Kl — > e/i is considered. In particular, we find the inequality 



^~4.xlO^ (5.5) 

If we assume the lowest M/^j-boson mass allowed from K^K^ mixing [^^], i.e., Mr = 1.5 
TeV, this then leads to mu < 45 GeV. This bound becomes much weaker for heavy Wr 
bosons, e.g., mr) < 0.5 TeV for Mr = 5 TeV. 

In the manifest LRSM, the largest effect in the decay amplitude due to the box graphs 
comes from the top quark, whereas charm-quark effects are in general sub-dominant. The 
relative behaviour of the charm- to top-quark contribution is given by the ratio 



mrV:,Vr.. 4 8 

' c/t 



^box ^ ""cycd^'CS - ^ ^ C5 Q) 



From Eq. ( [5. 61 ), we see that top-quark mass effects dominate over the charm-quark effects in 
the box diagrams, when mjv ^ Mr. The situation is very different for the K^K^ mixing, in 
which charm-quark yields the most dominant contribution in the loop |l^ . The reason for 
this crucial difference is that the one-loop box functions show up an enhanced behaviour, 
when the largest fermionic mass in the loop is much bigger than Mr. Obviously, in K^K^ 
mixing, the highest fermionic mass scale is set up by the top quark, which is much smaller 
than the V^/^-boson mass, whereas the heavy neutrinos can be much heavier than Wr in 
the decay amplitude for Kl e/i. In the manifest LRSM, we can estimate the relative 
contribution of the one-loop box functions for the two cases mentioned above as follows: 

Aox(i^° - K^) ^ V:aVcs , . 
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This discussion shows how heavy Majorana neutrinos with masses ttin ^ Mr are phe- 
nomenologically very significant for enhancing the decay rate of e/i. 

The decay Kl — > e/x can also proceed via a tree-level FCNC Higgs exchange. There- 
fore, it is most important to compare the FCNC self-energy contribution with that of the 
tree-level amplitude. In particular, we find the ratio 

~ « ~ qgO _ 380) x — | . (5.8) 



S(3(7 + 3/i) a^O{\)Ml ^ ' Ml 

It is easy to see that B{?)g + 3/i) is getting larger than iStree for M^ ~ ISM/j. On the 
other hand, the constraint that theory should remain perturbative at high energies gives 
rise to the unitarity bound in Eq. ( |3.10| ). As opposed to K^K^ mixing, the tree- level Higgs 
contributions to the decay Kl — > e/x are always comparable to the one-loop effects, and 
hence both must be included in the analysis. This constitutes another non-trivial feature 
for the decay Kl — > efi, which is to be examined carefully in the numerical estimates. 

In the numerical estimates, we consider two representative variants of left-right mod- 
els: (i) the manifest or pseudo-manifest LRSM {i.e. = = V), in which the Dirac 
mass matrix is still different from the charged- lepton mass matrix Mi, and (ii) the 
non-manifest LRSM, where = V but ^ V. Implementing the results of Ref. |22 



we choose natural relations and mass hierarchies among the kinematic parameters of the 
left-right models. In particular, we set all heavy Higgs particles to be degenerate, i.e., 
M^r = M^i = M^ and Mh = M^ {cf. Eq. ( p.6|) ). Already a small mass difference between 
the FCNC Higgs scalars would be sufficient to induce a large negative radiative shift in 
Rf, = T{Z — s> })b)/T{Z hadrons), which would be incompatible with the existing tight 
constraints on this LEP observable W\ 



To reduce even further the number of the many free parameters, we fix (s^")^ = 10 ^, 
^zT)^ = and the angle = 7r/4 in Eq. ( p.20| ). This choice of heavy-light neutrino mixings 



is in compliance with other low-energy constraints derived in Eqs. ( p.24|) and (|2.25|) . For 



simplicity, we also assume that the two heavy neutrinos, A^^i and N2, have the same mass, 
i.e., rriNi = = ^n, and the heavy neutrino mass uin- is varied from 100 GeV up 
to the perturbative unitarity limit 50 TeV which is compatible with the limit of Eq. 



(^). The input values for the entries of the CKM matrix V = are ||]: V^d = 0.224, 
Vcs = 0.975, Vtd = 0.015 and Vts = 0.048. In the non-manifest LRSM, the mixing matrix 
for the right-handed quarks is in general arbitrary and is only subject to phenomenological 
constraints. Here, we set = V^f = 1 and = V^f = 0, which leads to an enhanced 
charm-quark contribution to the decay amplitude, and in this respect, this model differs 
from the manifest LRSM. 
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Since we are interested in investigating the individual tree and one-loop contributions 
to the decay amplitude T{K^ e~^fi~), our numerical estimates for B{Kl e/i) are 
computed from the four different squared matrix elements: 

I. iT^reeP (solid), 

(5.9) 



II. \Ttrccr + 2^e{Tt,,ai:ioop) (dashed), 



III. ITi.ioopP (dotted), 

IV. |7^rcc + ^i-ioopP (dash-dotted). 

In Eq. ( [5. 91 ), the type of lines used in the plots for the four different predictions is specified 
within the parentheses. The results obtained from the expressions I and II are those that 
are based on a consistent loop expansion of the squared decay amplitude |T(A'° — * e+/i~)p. 
For comparison, we also give numerical estimates derived from the formulae III and IV. 

In Fig. 7, we display the dependence of B{Kl efi) on the heavy Majorana mass 
rriN for the selected values of the IVfj-boson mass Mfj = 1, 2, 5 and 10 TeV in the manifest 
LRSM. In agreement with Eq. ( |5.3| ), we find a strong dependence of B{Kl —>■ efi) on m^. 
This non-decoupling behaviour of the heavy neutrinos enters via the Dirac mass terms 
rriD of the see-saw type mass matrix in Eq. ( ^.9| ) JT^. We see that the pure one-loop 
contributions to the squared decay amplitude (dotted line) are always comparable to the 
respective tree- level results (solid line). In the manifest LRSM, the V^R-boson mass is 
severely constrained by the experimental limit on the K^K^ mixing. After including QCD 
corrections, the authors in find M/j ~ 2.5 TeV in the limit of very large FCNC Higgs 
masses. In all plots in Fig. 7, we take = 15Mfi, for definiteness. As can be seen from 
Fig. 7(c), for example, for sufficiently heavy Majorana neutrinos, B{Kl —>■ efi) may exceed 
the present experimental bound, which is indicated by an horizontal dotted line in all plots. 
Thus, for Mji ^ 2.5 TeV, we find that mjv ~ 10 TeV, for the set of the input parameters 
mentioned above. 

Fig. 8 gives the numerical predictions obtained in the non-manifest LRSM. We use 
the same input parameters as those of Fig. 7, except of the fact that we now consider the 
afore-mentioned matrix, which significantly suppresses the top-quark contribution in 
the loop. In this scenario, the chirality mass enhancement due to the heavy top quark is 
no longer applicable and pure loop-effects (dotted line) on e/i are at most 10% in 

comparison to the tree- level results for the experimentally accessible region. In fact, the 
tree-level effects on Kl — > e/i are more significant than the manifest LRSM case. The latter 
leads to tighter constraints on the heavy Majorana mass mjv. It is interesting to note that 
the non-manifest LRSM at hand does not invalidate the experimental bound due to K^K^ 
mixing. As a result, even relatively light heavy Majorana neutrinos with masses of few 
hundreds of GeV can account for possible lepton- flavour violating decays of Kl. 
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In Figs. 9 and 10, we examine the dependence of B{Kl efi) as a function of Mr in 
the manifest and non-manifest LRSM, respectively. We vary the V^R-boson mass from 0.5 
to 20 TeV, for different values of the heavy neutrino mass rriN- The FCNC Higgs scalar 
mass Mrp is also varied in linear dependence of M/j, as is explicitly shown in the plots. As 
is expected, we observe the generic feature of decoupling of a very large Wr mass from the 
theory. In the manifest LRSM (see Fig. 9), B{Kl efi) is difficult to measure for M/j > 7 
TeV, even if rather heavy Majorana neutrinos and FCNC Higgs bosons are assumed. In the 
case of non-manifest LRSM, the respective mass bound is slightly weaker, i.e., Mr > 10 
TeV, as can be seen from Fig. 10. 

It is now interesting to investigate the impact of both the FCNC Higgs-scalar and 
Wfl-boson masses on B{Kl —>■ e/i). To this end, we present exclusion plots in Figs. 11 and 
12 for the manifest and non-manifest LRSM, respectively. The allowed areas in the M^-Mr 
plane are determined from the inequality B{Kl e/i) < B^xp, where B{Kl — > e/x) has 
been calculated using the squared matrix elements in Eq. (|5.9| ) for different values of the 
heavy neutrino mass. The solid, dotted and dashed lines exclude the areas, in which the 
'm^v' label is not contained. Clearly, the tree- level results (solid line) are independent of the 
Wr mass, and are in qualitative agreement with Eq. ( p.9|) . In contrast, the results derived 
from the one-loop expression III in Eq. ( |5.9|) (dotted line) are insensitive to variations of 
M$. We include the constraints obtained from the squared matrix element III, since they 
will become significant beyond one loop. The dashed lines represent the limits obtained 
from the expression II in Eq. (|5.9|) , whereas results based on the squared matrix element 
IV are not shown. To a good approximation, the latter may be represented equally well 
by the combined graphical effect of all the three lines mentioned above. From Figs. 11 and 
12, we can see that a large region of the parameter space of the two typical LRSM's can be 
excluded by the experimental limit on the decay e/i. For the values of the light-heavy 

neutrino mixings (s^")^ = 10""^, {s'^lY = 0, a relatively low Wr mass, i.e., Mr < 1 TeV, 
cannot naturally accommodate the experimental data from Kl —>■ e/i and K^K^ mixing, 
without recourse to excessive fine tuning in the parameter space of the left-right models. 



6 Conclusions 

We have studied the lepton-fiavour- violating decay Kl cfi within the framework of 
SU(2)r(8)SU(2)/^®U(1)(b„l) models. Experiments searching for this decay mode are very 
important, since they can offer complementary limits on the parameter space of the left- 
right models together with limits obtained from the absence of fi ^ cy, ^ — e conversion 
in nuclei and fi eee, and from the observed K^Ks mass difference. 
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We have seen that one-loop box diagrams are not sufficient to provide a gauge- 
invariant result. Working in the gauges and adopting the well-defined OSS renor- 
malization, we have shown explicitly how gauge independence is restored in the transition 
amplitude, T{K^ —>■ e~^fi~), at one loop, when the corresponding self-energy and vertex 
graphs involving the FCNC Higgs scalars are taken into account. In fact, in the Feynman-'t 
Hooft gauge, the gauge-dependent complements of the box graphs may become dominant 
for a large range of the parameters. This novel aspect has not been addressed in detail in 
the existing literature before. 

Using gauge- independent analytic expressions for the one- loop matrix element, we 
have performed a systematic analysis studying the explicit dependence of B{Ki efi) on 
the various kinematic parameters. If we keep the enhanced light-heavy neutrino mixing 
s'2 ~ m£,/mM fixed and increase the high Dirac mass terms in the left-right model, which 
arise from the spontaneous breakdown of the SU(2) l gauge symmetry, we obtain a quadratic 
dependence of the one-loop transition amplitude on the heavy neutrino mass. In contrast, 
if we hold fixed below the perturbative unitarity bound but vary the Majorana mass 
scale mM, then the transition amplitude arising from the left-handed currents, e.g., from 
two bosons, vanishes in the limit niM oo, namely. All — > 0. As is expected on 
theoretical grounds, heavy neutrinos being singlets under SU(2)i will eventually decouple 
in this limit. Of course, these are known facts that have already been observed in [|1^ 
and subsequently discussed in many articles [|l^, p!5|JT6|JT7|l . Finally, our numerical estimates 
confirm the non-decoupling behaviour of the heavy neutrinos due to high SU{2)l Dirac 
mass terms. 

In LRSM's with enhanced light-heavy neutrino mixing, however, there exists an ad- 
ditional heavy-neutrino enhancement due to the simultaneous presence of left- and right- 
handed currents in the loop |ffOl . These chirality- changing currents give rise to the reduced 



amplitude B in the decay amplitude (3J-). Since tum is not a iso-singlet mass term under 
SU(2)ij, we have found that B tends to the non- vanishing expression iritmL)/ M'^ in Eq. 
( p.3|) for niM — ^ oo and fixed VEV vr. Clearly, mulvLi cannot be arbitrarily large but 
should satisfy the perturbative bound (|2.27| ). It may be worth emphasizing again that 
the observed non- decoupling behaviour of heavy neutrinos in the reduced amplitude i3 is a 
novelty in the left-right models and has no analogue in the SM with right-handed neutrinos 
discussed above. In fact, we have found that the new non-decoupling phenomenon can 
only originate from amplitudes, which violate chirality in the loop, thereby allowing for an 
active interplay between the left-handed mass scale and the right-handed one mu- 

Numerical estimates have shown that a significant part of the parameter space of the 
left-right models may be constrained due to the dominant contribution of B and other low- 
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energy data. Since these constraints depend in general on many independent parameters, 
we have systematized the presentation of our results in a number of plots in Figs. 7 - 12. In 
the numerical estimates we have studied the manifest as well as a typical non-manifest left- 
right model. To exemplify further the significance of such an analysis, we note in passing 
that for m^v ~ 10 TeV, s"]^ 10~^ and very heavy FCNC Higgs scalars, exclusion plots lead 
to Wr boson that must be heavier than about 3 TeV in the manifest left-right symmetric 
models (c/. Fig. 11(d)). 

In conclusion, if experiments at DA$NE or other future machines could substanti- 
ate a non-vanishing B{Kl — > eji) at the level of 10"^^ — 10~^^, this would be naturally 
accounted for within left-right models through the novel non-decoupling behaviour of the 
heavy neutrinos in the chirality-changing amplitudes. In the manifest left-right symmetric 
models, such an explanation would require relatively light VFij-boson masses with Mr 2 
TeV and heavy neutrinos of several TeV, for light-heavy neutrino mixings ^ 10~^, 
constraints which are in compliance with all other low-energy data. 

Acknowledgements. K.S. would like to thank Ben Grinstein and the UCSD physics 
department for their hospitality. The work of K.S. is supported in part by the Volkswagen- 
Stiftung. 
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A Feynman rules in the LRSM 



In this appendix, we list all the relevant Feynman rules obtained within the left-right 
models, which govern the interactions of the gauge and Higgs bosons with leptons and 
quarks, as well as the trilinear self-couplings of the bosons. In particular, we assume that 
the SU(2)i weak coupling constant is equal to the corresponding SU(2)r one qr, i.e., 
9l = 9r- 

With the assumptions and simplifications mentioned above, the trilinear couplings 
among gauge, Higgs, and would-be Goldstone bosons, which are only relevant for Kl — > e/x, 
are given by (all momenta flow into the vertex) 
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where sp is defined in Eq. ( ^.51 ) and are the FCNC Higgs bosons. 

The corresponding couplings of the gauge, Higgs, and would-be Goldstone bosons to 
the charged leptons and neutrinos can be read off from the Lagrangians: 
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where the mixing matrices and are defined in Section 2, and summation over 
repeated indices is imphed. 

Similarly, the interactions of the gauge, Higgs and would-be Goldstone bosons with 
the quarks may be obtained from the Lagrangians 
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where and are the corresponding left and right mixing matrices in the quark sector. 



B One-loop analytic results 



Here, we will present the analytic expressions of all one-loop results. To this end, we first 
define the following loop functions: 
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Equipped with the above loop functions and the Feynman rules given in Appendix 
A, we can calculate the individual contribution of the box diagrams, shown in Figs. 2-5, 
to the decay amplitude of Kl —>■ e/j,. More explicitly, we have 
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where p = M^/M^, A^ = M^/M^, and A^ and A^ are defined after Eq. ( ^) . 

As has been shown in Section 4 however, the sum of the box contributions to B 
is not gauge independent by itself. One has to include vertex and self-energy graphs, 
which are depicted in Figs. 3 and 4. These graphs will restore gauge invariance, after 
renormalization has been taken into account (see Section 4). Since we are interested in 
finding the corresponding gauge-dependent complement part of the box contributions to 
B in the Feynman-'t Hooft gauge, we keep only those terms that do not vanish in the 
formal mass limit M^r,i — > oo. In Section 4, we have also seen that there will be an implicit 
dependence of the loop integrals on the FCNC scalar masses, which enters via the OSS 
renormalization. Furthermore, we take both the FCNC scalars to be exactly degenerate, 
i.e., M^r = M^i^ = M^, which is a good approximation (see also Eq. ( |2.(j| )). In this way, it 
is not difficult to find from the would-be Goldstone-boson graphs 
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Since A/i — = 0{1), as has been discussed in Section 2, the reduced amphtudes B^ic+M), 
B[Ae + 4/) and B{Ag + Ah) behave as l/A^ and 1/A^. Therefore, these contributions may 
safely be neglected in the numerical estimates. 

The OSS renormalization scheme adopted in Section 4 gives rise to the subtracted 
loop functions Cq^ and which are defined in the following way: 

5^2(M^m?,m^) = So(0,m?,m^)-fio(M^m?,m^) + M2s;(M^m?,m^),(B.24) 
Co^^(0,M^O,m?,m2,m^) = Co(0, 0, 0, m?, m^, m^) -Co(0, 0, m^, m^, m^) . (B.25) 

The loop functions Bo{M'^ , mf, m"^) and Co(0, M^, 0, m^, ml, mg) are the usual Passarino- 
Veltman functions defined in Ref. Their explicit form is given by 
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where A(x, 2;) = (x — y — z)'^ —^yz, cosh "^z = ln(z + V — 1) and Cuv is an UV constant 
which drops out after renormalization. Furthermore, in Eq. ( p.27|) , we have defined 
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and the dilogarithmic function 
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Figures 7-12 





Fig. 7: B{Kl efi) versus mjv for different M/j values in the manifest LRSM 
(the meaning of the various Hues is given in the text). 
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Fig. 8: B{Kl — e/i) versus mjv for different Mr values in the non-manifest LRSM 
(the meaning of the various hnes is given in the text). 
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Fig. 9: B{Kl — e/i) versus Mr for different mjv values in the manifest LRSM 
(the meaning of the various hnes is given in the text). 
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Fig. 10: B{Kl — e/i) versus Mr for different mjv values in the non-manifest LRSM 
(the meaning of the various hnes is given in the text). 
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Fig. 11: Exclusion plots for the parameters M^, Mr, obtained from the constraint 
B{Kl e/i) < 3.3 X 10^^^ in the manifest LRSM (the meaning of the 
various lines is given in the text). 
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Fig. 12: Exclusion plots for the parameters M^, Mr, obtained from the constraint 
B{Kl e/i) < 3.3 X 10^^^ in the non-manifest LRSM (the meaning of the 
various lines is given in the text). 
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